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Josephson junction arrays can be used as quantum channels to transfer quantum information 
between distant sites. In this work we discuss simple protocols to realize state transfer with high 
fidelity. The channels do not require complicate gating but use the natural dynamics of a properly 
designed array. We investigate the influence of static disorder both in the Josephson energies and 
in the coupling to the background gate charges, as well as the effect of dynamical noise. We also 
analyze the readout process, and its backaction on the state transfer. 



The transmission of a quantum state through a channel 
between distant parties is an important issue in quantum 
communication. In optical systems photons can be trans- 
ferred coherently over large distances jjj. However, it is 
also highly desirable to have similar protocols for quan- 
tum information transfer in solid-state environments. A 
possible solution would be to interface solid-state quan- 
tum hardware to optical systems 0] . Another possibility 
is to use flying qubits, i.e. to transfer the physical qubits 
along leads U • Inspired by the paper of Bose 0] the idea 
of our work is to construct a genuine quantum transmis- 
sion line using a Josephson junction array. 

Recently, a spin chain with ferromagnetic Heisenberg 
interactions has been proposed for quantum communica- 
tion 4] . It was shown that Heisenberg chains can be used 
to transfer unknown quantum states over appreciable dis- 
tances 10^ lattice sites) with high fidelity Hllli,[3|. 
By preparing the state to be transferred at one end of 
the chain and waiting for a well-defined time interval, one 
can reconstruct the state at the other end of the chain. 
Even perfect transfer could be achieved over arbitrary 
distances in spin chains 0. Quantum state transport 
through harmonic chains was considered in Ref . 0. 

Josephson qubits are among the most promising candi- 
dates as building blocks of quantum information proces- 
sors In this Letter, we extend their application 
range to quantum communication and show that a one- 
dimensional Josephson array is a natural transmission 
line for systems with superconducting charge qubits. We 
calculate the transmission fidelity and investigate the ef- 
fect of static inhomogeneities and dynamical noise. We 
also analyze the readout process by a single-electron tran- 
sistor (SET) at the end of the array. 

The model that we want to study is schematically il- 
lustrated in Fig. n and described by the Hamiltonian 



H = H}3 + Hqp + Hcoup, where 
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Hjj = ^ ^(Q,; - Qxi)C.^j ^ {Q-j - Qxj) cos (f>,.,+i 

ij i 

(1) 



FIG. 1: Dashed box: one-dimensional Josephson array pro- 
posed for the transmission of quantum states. The crossed 
rectangles denote the Josephson junctions between the is- 
lands. The state prepared on the left-most island is transfered 
to the right-most island by the time evolution generated by 
the Hamiltonian. Left part: Cooper-pair box (charge qubit) 
used to prepare the state. Right part: SET transistor used as 
measurement device. 



is the Hamiltonian of a one-dimensional Josephson junc- 
tion array 12] of length L, and (pi^i+i = (pi — 4>i+i- The 
other terms of the Hamiltonian describe the measurement 
apparatus and will be discussed later. The charge Qi 
and phase pi are canonically conjugated. The first term 
in Eq. |^ is the charging energy, Cij is the capacitance 
matrix; the second is due to Josephson tunneling. An ex- 
ternal gate voltage Vxi gives a contribution to the energy 
via the induced charges Qxi = 'ieqxi = VxiCa. This ex- 
ternal voltage can be either applied to the ground plane 
or unintentionally caused by trapped charges in the sub- 
strate (in this case Qxi will be a random variable). We 
assume that each island is coupled to its nearest neigh- 
bors by junction capacitances C and to the ground by a 
capacitance Cq. In this case, the charging interaction has 
a range given by C /Cq in units of the lattice spacing 
of the array In the following we put h= ks = 1- 

In the charge regime c^C^q :s> Ej, the system is 
approximately described by only two charge states for 
each island. The chain Hamiltonian i/jj is equivalent to 
an anisotropic XXZ spin- 1/2 Heisenberg model [THl3|. 
the Josephson chain is thus different from the XY and 
Heisenberg cases 4]. It is characterized by a strong 
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anisotropy between the z-direction and the xy-plane. 
Moreover the z-coupUng has a range which depends on 
the electrostatic energy and can extend over several lat- 
tice constants. 

At i = 0, the chain is initialized in the state |'0o) = 
l^/i, 000...0}, where |0) (|2)) denotes the state of an is- 
land without (with) an excess Cooper pair, and — 
cos(6'/2)lO) + e''^sin(6'/2)|2) is the state that has been 
prepared in the left-most island. This initial state is not 
an eigenstate of the Hamiltonian, it will evolve as a func- 
tion of time. In fact, as the total charge Q = Yli Qi is 
a conserved quantity, the dynamics is restricted to the 
L + 1-dimensional space H = 7io (B 'H2 of total charge 
zero, Hq, and charge two, H2 = span{[j)}, where 
1 < j < L is the state with an excess Cooper pair on the 
j-th site. In this basis the Hamiltonian reads 



E, 



((l-<5,x)|j + l) + (l-<5,i)|j-l)) . (2) 



We first calculate the fidelity of transmission and the 
time required for the transfer of information as a function 
of the coupling constants of the Josephson chain. The 
quality of the transmission is quantified by the fidelity of 
the (mixed) state pl of the right-most island (site L) to 
the initial state 
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This definition gives the fidelity averaged over all possible 
initial states on the Bloch sphere, 1/2 < Fl < I. 

The fidelity is a strongly oscillating function of time. 
Only at well-defined times the state is transferred faith- 
fully through the chain. This does not necessarily corre- 
spond to the time in which a Cooper pair has been trans- 
ferred, since also the relative phases of the state have to 
be reconstructed. In Fig. [21 we show the value of the first 
fidelity maximum and the time at which it is reached as 
a function of the length L of the array and for different 
values of the ratio C/Cq- For the parameters considered, 
the fidelity is never smaller than 75%. For longer chains, 
or if the condition Co ^ C is released, the first maximum 
of the fidelity is considerably reduced. Another option is 
to fix a threshold for the fidelity of transmission and seek 
for the first local maximum above the threshold. The 
time at which these maxima occur increases exponen- 
tially with the chain length. The value of the fidelity 
does not necessarily decrease on increasing L, and for 
larger arrays a higher fidelity can be achieved (although 
at larger times), see Fig. 13 The results of Figs. |2 Elare 
encouraging since they indicate that faithful state trans- 
mission using Josephson chains is already possible with 
present-day technology. 



Since experimental arrays are never completely ho- 
mogeneous, we now consider the case in which a small 
amount of static disorder is present. In general, imperfec- 
tions will reduce the fidelity. In Fig. 0|we show both the 
effect of bond disorder (Josephson coupfings distributed 
around an average value) and site disorder (mimicking 
the effect of static background charges and/or different 
capacitances). The effect of charge disorder appears to 
be more disruptive: this is because additional frequencies 
enter the dynamical evolution making the reconstruction 
of the additional wave- function more difficult. 

Dynamical fluctuations play a different role. They 
arise from gate-voltage fluctuations and are described 
by adding stochastic terms to the gate voltages, qxi — > 
Qxi + £,i (t) in the Hamiltonian in Eq. . Here we choose 
a very simple model and assume the (t) to be indepen- 
dently gaussian distributed: {£,i{t)) = 0, {£,i{t)^j{t')) = 
^dijS{t — t'). Nevertheless, due to capacitive couphng 
between separated sites, such stochastic factors result in 
correlated stochastic terms in the effective Hamiltonian 
Eq. ©, i?noisc|j) = Hjj\j) - 2Ej{t)\j), where the zero- 
averaging gaussian functions (t) , are uniquely fixed by 
{E,{t)Ej{t')) = j[{C-^)%S{t - t'). Averaging out the 
stochastic terms leads to the master equation for the den- 
sity matrix of the chain in the space Ti, 



P = -i[H.]j,p] 



(4) 



(QiQjP - 2QipQj + pQiQj) 



where the operators Qi projected on the space Ti. are 
Qi = 2e|i)(i|. The state of the system develops into 
a completely incoherent mixture in which any charged 
state is equally probable, pii{t 00) — {1/L)pii{t = 0). 
The averaged fidelity as defined in Eq. Q is reduced to 
Foo = 1/2-1- 1/(6L), corresponding to an almost unfaith- 
ful transmission. The time dependence of the fidelity in 
the noisy system is presented in Fig. O where it is com- 
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FIG. 2: Maximum value of the fidelity as a function of the 
length of the chain for two different values of C/Cq <^ 1 and 
(2e)^/(£',/Co) — 10. Inset: time at which the maximum is 
reached. 
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pared with the fidehty in the absence of noise. The peaks 
of the fidehty are not smeared out by noise. The dom- 
inant effect of the coupling to the environment is the 
relaxation of the fidelity amplitude towards the station- 
ary value (independent on the initial state). Numerically, 
such relaxation takes place on a characteristic time scale 
^ 1/(^7). Thus, to observe high values of the fidelity it 
is important to have a maximum at a short time. As a 
consequence, C/Cq ^ 1 is preferable, and the condition 
7 ^ Ej 1 1? is required to have a high value for the first 
maximum of the fidelity. 

Finally we discuss how the fidelity can be measured 
in a practical setup. To do this, we assume that the 
right-most island (site L) is part of a SET transistor. 
We therefore specify the effective coupling Hamiltonian 
between the right-most island and the leads 0|, 



b—u,d.L k.(7 



(5) 



b—u,d.L 

- ^ Jb cos{(f>L - 'fib) - ^ VtQb , (6) 

h—u,d h—u,d 



where 7, 7^ are annihilation and creation operators of 
quasiparticles in the grain L and in the leads, u and d (see 

Fig- HI- The operator Xb = Y.k,q,a Tqkll^-hlq^L describes 
quasiparticles tunneling into the grain with an associated 
charge increasing e^'^'^^^'''''^/^. is the total charge 
entering the chain from the up or down reservoirs. We 
assume non- vanishing quasiparticle tunneling only across 
the upper junction {eVd = 0, eT4 = —eV « — 2A) and 
conversely we allow coherent Cooper pairs tunneling only 
in the lower junction {Ju = 0, Jd = J Ej). Gate 
voltages are chosen so that the SET is off resonance and 
we can therefore neglect the stationary current through 
the SET due to the Cooper-pair quasiparticle cycle [I^ . 
The measurement device modifies the dynamics of 
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FIG. 3: Maximum value of the fidelity as a function of the 
length of the chain for three different values of C /Co ^ 1 and 
{2e)^ /{EjCo) ~ 10. Inset: time at which the maximum is 
reached. 
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FIG. 4: Fidelity as a function of time for an array of length 
L — 7, {2e)^ /{EjCo) = 10 and C = 0, i.e. a junction capac- 
itance much smaller than the ground capacitance. Disorder 
parameters: relative variance AEj/Ej — 0.1 for bond disor- 
der, absolute variance AQx/2e — 0.025 for site disorder. 



Cooper pairs on the chain and requires taking into ac- 
count quasiparticle excitations on the L-th site of the 
chain. By neglecting quasiparticle tunneling we would 
have a coherent dynamics for the charges in the chain 
described by the Hamiltonian Hq — H{Tq^. 0). Trac- 
ing out the quasiparticle degrees of freedom results, in- 
stead, in an incoherent dynamics described by a mas- 
ter equation for the reduced density matrix p of charges 
in the chain In the basis of eigenstates of Hq, 

Ho\a) — Ea\a)V?\, the master equation reads 



Pap{t) = -i{a\ [Hq,p] 1/3) 



RaPiivPiiv ■ (7) 



The prime indicates that the sum has to be performed 
over states with energies such that | Ea — E/s — Ef^ + E^\ ^ 
1/Ai, At being the time over which the coarse-graining 
implicit in Eq. (0 takes place. 

As we are interested in the time evolution over short 
times, let us discuss in some detail our approxima- 
tions. We first assume that J <C Ej, so that, in 
evaluating the kernel R, we neglect the Josephson cou- 
pling to the leads J^l. In this case the spectrum of 




FIG. 5: Fidelity versus time in presence of gate voltage fluc- 
tuations (full line). Dashed line: noiseless case. L — 7, 
e^/iEjCo) = 10, C/Co = 0.1, 7 = 0.01. 
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FIG. 6: Full line: time dependence of the current (in units of 
e/T) through the SET. Dashed line: fidelity of the isolated 
chain. F = 0.05, all other parameters as in Fig. |5] 

is {£'0, i?i/2, i?j\7}, M = 1,...,L. Due to the en- 
ergy scale separation Ej^^ — Eff ^ Ej <^ 1/At < 
Em - El Em - Eq <t: eV, the sum in Eq. Q 
mixes population and coherences of the density ma- 
trix only in the subspace of M states. In this case 
the coarse-grained dynamics of Eq. (TJ can resolve the 
time scales of order l/Ej that we are interested in. 
In this approximation, the only non-vanishing terms of 

the kernel R are found to be i?nnii — — ^^iiii — 

^^22 2222 

{Rmnpq} ^ ^fef.pAjvQ + Sn^qXmp], where Xab = 
T{A\L){L\B). In these expressions we approximated 
r ~ dsexp{ieVs){Xu{s)Xl{0)) ~ dsexp(i(el^± 
(Em — Ei)s){Xu{s)Xl{0)) as a consequence of the sep- 
aration of energy scales discussed above. We also ne- 
glected all other exponentially small (~ e"*^^/"^) rates. 

Finally let us address the proposed measurement pro- 
tocol. It consists in disconnecting the right-most site 
from the rest of the chain at a time t* ^ 1/F and in 
measuring the time-integrated current through the SET 
I=Tir ^(^)'^^ = f It^ i{T)dT+0{Tt*) where T is the 
time between two pulses: it is the largest time scale in 
the system. The instantaneous particle current is I{t) — 
^{pLhit) + /5i 1 (t)). The last term 0{Tt*) in the current 

corresponds to /q i{T)dT = {pLL{t) + Pkk{t)) <Tt* . 
As at time t > t* the SET is disconnected from the rest of 
the chain and is out of resonance, the measured current 
is '-^ ^{2pLL{t*) + pi This measurement scheme 

does not provide a tomography for the state of the right- 
most site: the measured current does not depend on the 
coherences of pL{t*), to which the fidelity is sensitive. 
Nevertheless, the peaks in the current correspond exactly 
to the maxima of the fidelity as shown in Fig.|B| The cur- 
rent decay in time due to quasiparticle tunneling happens 
on a time scale ^ 1/F irrespective of the length of the 
chain. Therefore, our measurement scheme can be used 
also for long chains, the main constraint are disorder and 
gate voltage fluctuations. In this sense the current mea- 
surement allows to check the theoretical prediction for 



the fidelity of state transfer. 

In conclusion, we have proposed to use a Josephson 
junction chain as a solid-state quantum communication 
channel. We have analyzed the read-out process at the 
end of the channel and shown that the fidelity can be 
directly measured in a SET device. We have also consid- 
ered the influence of static disorder and dynamical noise. 
Present-day technology should allow the realization of 
quantum channels of the type described here. 
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sion in a setup using SQUID loops coupled to resonators. 
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